Topographic and electronic contrast of the graphene moire on Ir(lll) probed by 
scanning tunneling microscopy and non-contact atomic force microscopy 
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Epitaxial graphene grown on transition metal surfaces typically exhibits a moire pattern due to 
the lattice mismatch between graphene and the underlying metal surface. We use both scanning 
tunneling microscopy (STM) and atomic force microscopy (AFM) to probe the electronic and to- 
pographic contrast of the graphene moire on the Ir(lll) surface. STM topography is influenced by 
the local density of states close to the Fermi energy and the local tunneling barrier height. Based 
on our AFM experiments, we observe a moire corrugation of 35±10 pm, where the graphene-Ir(lll) 
distance is the smallest in the areas where the graphene honeycomb is atop the underlying iridium 
atoms and larger on the fee or hep threefold hollow sites. 

PACS numbers: 68.37.-d 68.37.Ps 68.65.Pq 



Epitaxial graphene can be grown on many transi- 
tion metal surfaces using chemical vapor deposition 
(CVD)i~— This process gives ready access to high- 
quality, large scale, graphene monolayers on surfaces 
where graphene growth is self-terminating (e.g. Cu, Ir, 
Pt)J-2ai These layers can be characterized by surface sci- 
ence techniques and if necessary, transferred onto other 
substrates for further processing. The different metal 
surfaces can be coarsely classified based on how strongly 
the graphene layer interacts with the underlying metal 
substrate^ For example, Ir(lll) and Pt(lll) surfaces in- 
teract weakly with the graphene layer and consequently, 
graphene still exhibits linear Dirac-like dispersion char- 
acteristic of isolated graphene.—"— On the other hand, on 
Ru(0001) and Ni(lll) surfaces, the graphene band struc- 
ture is strongly modified. 5,8 While the CVD growth oc- 
curs epitaxially, the lattice mismatch between graphene 
and the metal substrate gives rise to a moire pattern that 
is observed on most metal surfaces (notably Ir(lll)^"— , 
Rh(lll) 13 , Ru(OOOl)^--^ and Cu(lll)H). 

This moire pattern can be readily observed by scan- 
ning tunneling microscopy (STM) ^ 14 i 17 However, STM 
images do not directly probe the topography of the sur- 
face; instead, the STM tip traces constant integrated lo- 
cal density of states (LDOS) surfaces at energies close 
to the Fermi level . 18 i 19 This causes the contrast and ap- 
parent corrugation of the graphene moire on Ir(lll) to 
depend on the STM imaging conditions. 9,10 It is not a 
priori clear which STM images correspond to the actual 
topography of the surface. 

The use of a quartz tuning-fork force sensor in the 
QPlus configuration has made it possible to carry out 
non-contact atomic force microscopy (nc-AFM) in the 
frequency modulation mode with small tip oscillation 
amplitudes. This allows concurrent STM experiments, 
where the performance of the STM mode is not compro- 
mised by the tip oscillation or small force constant of the 



AFM cantilever i 18 i 2Q ~ — We have used this technique and 
performed both low-temperature AFM and STM mea- 
surements on epitaxial graphene monolayers on Ir(lll) 
aimed at understanding the contributions of actual to- 
pography, charge transfer giving rise to local variations in 
the tunneling barrier height and contact potential differ- 
ence, and variations of the LDOS on the observed moire 
pattern. These techniques give independent information 
on the surface topography, which allows separating elec- 
tronic and topographic effects. 

The graphene was grown on Ir(lll) by CVD from 
ethylene.— The Ir(lll) surface was first cleaned by re- 
peated cycles of 3 kV Ar + -ion sputtering at room tem- 
perature followed by flashing to 1400 K and annealing at 
1200 K. After the last annealing cycle, the sample was 
first annealed for 3 min at 800 K in 1 • 10~ 7 mbar Oo, 
then flashed to 1400K before starting the CVD process. 
In order to ensure the formation of a full graphene mono- 
layer, the sample was exposed to 5 • 10 ~ 7 mbar ethylene 
at 1250 K for 100 s. The sample was then transferred ex 
situ to the low-temperature STM/ AFM system (LT-STM 
equipped with QPlus force sensor, Omicron Nanotechnol- 
ogy Gmbh). Prior to the STM/ AFM measurements, the 
sample was cleaned by heating to 800 K. All STM/ AFM 
experiments were carried out in ultra-high vacuum (base 
pressure < 10 -10 mbar) and at low-temperature (T = 4.7 
K). The QPlus sensor used for the frequency modulation 
nc-AFM experiments had a spring constant k of 1800 
N/m, resonance frequency /o of ca. 24 kHz and a quality 
factor of 18800. We used Ptlr tips and the tip oscilla- 
tion amplitude was set to 5 A. Bias voltage (Vbias) was 
applied on the sample with respect to the tip. d//dVbias 
and dl/dz signals were recorded with a lock-in ampli- 
fier by applying a small sinusoidal variation to the bias 
voltage or the z-piezo position, respectively. 

Figure QJi shows a constant-current STM topography 
image of a graphene monolayer on Ir(lll). In addition 
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FIG. 1. (Color online) (a) Constant-current STM topography image of epitaxial graphene on Ir(lll) acquired at Vbi as = 0.3 V 
and set-point current of 0.3 nA. The line indicates the moire unit cell and the three inequivalent areas within it are denoted 
by A, B and C. (b) High-bias STM image taken at 0.45 V / 1 nA showing the inverted moire pattern, (c) Constant-current 
d//dVbias and dl /dz maps recorded at a bias of 0.05 V. (d) Constant frequency shift nc-AFM image with Af = —45 Hz and 
Vbias = 0.01 V. (e) Average current over the tip oscillation cycle measured simultaneously under AFM feedback. 



to the atomically resolved hexagonal graphene structure, 
a moire pattern with a period of 2.5 nm is clearly visible. 
It has been shown previously that this superstructure 
preferentially orients along the atomic rows of graphene, 
which is also the case in Figure [TJi. 10 The apparent peak 
to peak (p-p) corrugation of the moire pattern is 50 pm. 
The unit cell of the moire is indicated by the solid line, 
and the three inequivalent areas (with respect to the reg- 
istry with the Ir(lll) lattice) are indicated by A, B, and 
C. It has been previously suggested that they correspond 
to areas where the graphene honeycomb is centered on 
the underlying Ir atoms (atop, A), or on the fee (B) or 
hep (C) threefold hollow sites We observed that the 
apparent corrugation depends on tip conditions and sys- 
tematically on the bias voltage, in agreement with earlier 
results. 10 STM images at a higher bias (> 0.5 V, Figure 
[TJd) exhibit inverted moire contrast (region A becomes 
bright) compared to low-bias images. 

In addition to standard STM imaging, we can get fur- 
ther information on the local electronic properties by 
mapping out the LDOS (oc dJ/dVbias) and the tunnel- 
ing decay constant k (oc dl/dz) signals in the constant- 
current mode. These quantities vary over the moire pat- 
tern as shown in Figure QJ. Both LDOS and n are lower 
at region A of the moire. 

As STM imaging is clearly influenced by electronic ef- 
fects, we carried out low-temperature nc-AFM experi- 
ments in the constant frequency shift mode to probe the 
surface topography of graphene monolayer on Ir(lll). 
Figure [T]i shows a typical AFM topography image (the 
image was not recorded on the precise location of the im- 
age shown in Figure [TJi) . In addition to the contrast on 
the atomic scale, we obtain a moire pattern with similar 
contrast as in STM imaging and an apparent corrugation 
of ca. 30 pm. Careful inspection of Figure [l]i reveals 
variations in the atomic scale contrast. On the bright 
areas of the moire pattern, the carbon atoms are im- 
aged as depressions, in line with the earlier atomically re- 
solved images of carbon nanotubes^ The observed con- 
trast changes on the dark areas, indicating that the tip- 
graphene distance is different on the different regions of 
the moire pattern (see below). 



The simultaneously measured tunneling current dur- 
ing the nc-AFM imaging is shown in Figure QJ. It again 
shows both atomic and moire contrast, where the low 
current regions are aligned with the depressions in the 
topographic image. Note that there is a shift between 
the AFM image and simultaneously measured tunneling 
current. It is likely that the tip has an impurity (atom) 
that does not contribute to the current but has an ef- 
fect on the measured frequency shift. Alternatively, an 
asymmetric tip apex can cause the shift between AFM 
and average current images. 2 - We find the same quali- 
tative moire contrast with different tips and on different 
locations of the sample. 

The tip-sample interaction causes a shift Af in the 
resonance frequency of the cantilever. At small tip os- 
cillation amplitudes, the measured detuning is directly 
proportional to force gradient, Af = —fo/(2k)(dF ts /dz), 
where F ts is the total interaction force between tip and 
sample. 18 Different forces contribute to F ts , the most rele- 
vant in our experiment being quantum mechanical forces 
between the tip apex and the surface (Pauli repulsion, 
chemical bonding) , vdW interactions between the tip and 
graphene and the tip and the Ir substrate, and electro- 
static forces j 18 i 19 i 23 AFM topography might also be af- 
fected by chemical inhomogeneity of the surface (different 
regions of the graphene moire are known to have different 
chemical reactivities ^ 13 i 14 ) . However, we have observed 
the same qualitative moire contrast with different tip ter- 
minations, consistent with the expected weak interaction 
between graphene and the Ir(lll) surface. 

Experimental results based on photoelectron spec- 
troscopy (ARPES and XPS) and ab initio calculations 
show that the interaction between graphene and iridium 
is weak££&2! Theoretical calculations give an average 
graphene-Ir(lll) distance of about 3.9 A (GGA-DFT) or 
3.42 A (LDA-DFT). 9,27 It is well-known that GGA un- 
derestimates and LDA overestimates binding in systems 
where vdW interactions are important. Recent DFT cal- 
culations using vdW-corrected functionals have found a 
binding distances of 3.6-3.7 A for graphene on weakly 
interacting metals (e.g. Pt). 28 Despite the large bind- 
ing distance, it is important to realize that the vdW 
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FIG. 2. (a) Schematic representation of the tip movement 
over graphene. (b) Variables used in Eq. [T] 



forces between the tip and the sample are sufficiently 
long range to include contributions from the iridium sub- 
strate. In the attractive regime, the background vdW 
from the Ir substrate results in increased attraction in 
the area A of the moire, which causes the AFM feed- 
back to increase tip-sample distance in order to keep Af 
constant. Hence, the AFM corrugation underestimates 
the real topographic corrugation of the graphene moire 
as illustrated in Figure [2ji. 

We will now model this effect within the small ampli- 
tude approximation (detuning proportional to force gra- 
dient). This approximation is valid if the force gradient is 
roughly constant throughout the oscillation cycle of the 
tip. We use relatively small oscillation amplitudes and 
consider here only the vdW interactions that are rele- 
vant for the AFM observation of the moire pattern on an 
otherwise chemically homogeneous surface. The chemi- 
cal interactions between the tip apex and the surface only 
act at very short distances (much smaller than our tip os- 
cillation amplitude) whereas the long-range electrostatic 
force remains almost constant over the oscillation cycle 
assuming that neither the tip nor the surface is charged. 
Consequently, these forces make only a minor contribu- 
tion to the observed AFM response. 

The total vdW force felt by the tip (modeled as a 
paraboloid z = x 2 /(2R), where R is the tip radius) 
can be calculated by integrating the vdW potential 
w v dw = — 4e(<j/r) 6 over the tip and Ir bulk, and tip 
and two-dimensional graphene layer. 19 We assume that 
the vdW interaction between the tip and Ir substrate is 
not screened by the graphene layer. Hence, we obtain 
an estimate of the upper limit of the background vdW 
contribution. The derivative of the total force is then 
proportional to the detuning Af of the tip 




dtip-gr^HOPG I 
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where Ai = 4:7r 2 €ip 2 af is the Hamaker constant and dtip-ir 
and dtip.gr are the tip-Ir and tip-graphene distances cor- 
responding to the midpoint of the tip oscillation cycle. 
In the case of a two-dimensional layer, the vdW force de- 
pends on the surface atom density rather than the volume 
density. We take this into account by using the HOPG 
Hamaker constant Ahopg and the layer density /hopg- 
In Eq. [I] the first term is the detuning caused by the 
tip-Ir vdW force and the second term the tip-graphene 
vdW force. The relation of these terms is illustrated in 
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FIG. 3. (Color online) Theoretical modeling of the vdW forces 
in the AFM tip-graphene-Ir substrate system, (a) The contri- 
butions of graphene and iridium substrate to the total Af. (b) 
Tip-graphene distance as a function of graphene-Ir distance 
calculated with Af = —45 Hz. The corresponding slope is 
plotted on the right ?/-axis. (c) Slope s of the dtip-ir vs. dti P - gr 
curve for different tip radii and different Af (indicated in the 
figure). Panel (a) and (b) are calculated with R = 10 nm. 
Ahopg = 3.42 • 1(T 19 J, Atip = Ai T = 4 • 1(T 19 J 29 ' 30 , and 
^hopg = 1/335.4 pm were used in all calculations. 



Figure [3^, where we plot them (and the total Af) as a 
function of the tip-graphene distance. The contribution 
from the Ir substrate increases and even becomes the 
dominant term at large distances. 

We have solved the tip-graphene distance dtip-gr from 
Eq. [1] numerically as a function of the graphene-Ir dis- 
tance (dtip-ir — <^tip-gr) (Figure ^jp) . dtip-gr changes almost 
linearly over a reasonable range of graphene-Ir distances. 
The slope s of this curve represents the underestimation 
of corrugation on the moire due to the background vdW 
forces from the Ir bulk. Thus the real geometric cor- 
rugation of the moire is given by a rea i = (1 — s)a afm- 
The Hamaker constants have a fairly small effect on the 
slope whereas R and Af have quite a significant effect. 
This is illustrated in Figure [3b, where we plot the av- 
erage slope in the range of d ZT .\ T between 3-4 A. This 
effect can be easily understood by the fact that Af and 
R determine the absolute tip-graphene distance, which 
governs the proportion of the vdW force from the bulk 
Ir with respect to the total force. While our estimation 
of dtip-gr depends on R, which is difficult to estimate in- 
dependently, we can also estimate dtip-gr based on the 
simultaneously recorded tunneling current. Taking into 
account the tip oscillation 18 and using a measured value 
of the tunneling decay constant k ~ 0.55 A -1 , we can ex- 
trapolate the distance to point-contact. This procedure 
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gives average tip-graphene distance of 7.2 A in Figure [T]i. 

The data shown in Figure [3h shows that even though 
the background vdW from the Ir bulk affects the ap- 
parent AFM corrugation, the effect is rather small. The 
correction factor is 10-20% for reasonable tip radii and at 
sufficiently negative A/, which yields 35±10 pm as our 
estimation for the actual moire corrugation. It should be 
noted that this model does not take the actual shape of 
the corrugation of the moire into account. It only cor- 
rects the measured corrugation by the background vdW 
from the Ir bulk. Due to the relatively equal scale of cur- 
vature of the tip and the moire, the corrugation may be 
slightly further underestimated by AFM measurements. 

The increase in the tip-graphene distance caused by 
the increased background vdW interaction between the 
tip and the Ir substrate in the area A of the moire reduces 
the tunneling current as seen in Figure However, if 
we use the estimated changes in dtip-gn we obtain a cur- 
rent variation of only ca. 10 %. On the other hand, the 
experimentally measured variation is much larger, ca. a 
factor of 2. Apart from the dtip-gr? the tunneling cur- 
rent is influenced by the LDOS and the decay constant 
k. Both of these quantities vary over the moire pattern 
as shown in Figure QJ. These quantities have opposite ef- 
fects on the tunneling current: the larger the LDOS, the 
larger the current. On the other hand, the larger the de- 
cay constant, the smaller the current as it is proportional 
to exp(— 2ftdtip-gr)- Our observations then imply that at 
small bias, the reduced LDOS in the region A of the moire 
is (mostly) responsible for the reduced tunneling current 
in the simultaneously measured tunneling current images 
under AFM feedback. This conclusion naturally does not 
hold for increased bias 0.5 V) where the STM contrast 
of the moire pattern is inverted. 

We can relate the STM to AFM results by switch- 
ing in situ back and forth between STM and AFM feed- 



back. Qualitatively, the moire contrast is the same be- 
tween AFM and low bias STM images (dark depressions 
in a bright background). We do not have a direct mea- 
sure of the registry between the moire unit cell and the 
underlying iridium lattice. However, comparison of our 
STM results with the STM-based graphene adsorption 
site determination- ! 10 ! 27 relates the regions of the moire 
unit cell to areas where the graphene honeycomb is cen- 
tered atop the underlying Ir atoms (A), or on the fee 
(B) or hep (C) threefold hollow sites. Thus, our STM 
and AFM measurements seem to imply that graphene- 
Ir(lll) distance is the smallest on atop sites (region A), 
and larger on fee and hep sites (regions B and C). This is 
surprising and in contrast to other graphene-metal sys- 
tems. 

In conclusion, we have carried out simultaneous low 
temperature AFM and STM experiments on an epitax- 
ial graphene monolayer on the Ir(lll) surface. These 
experiments shed light to the structure of the graphene 
moire on the Ir(lll) surface. While STM experiments 
are dominated by electronic effects, nc-AFM provides a 
qualitatively correct image of the surface topography. A 
more quantitative estimation of the moire corrugation 
based on the AFM experiments would require account- 
ing for the background vdW interaction between the tip 
and the metallic substrate. Although in the present case 
of graphene on Ir(lll) the background effect is small, it 
has to be considered in principle for quantitative topogra- 
phy of atomically thin two-dimensional layers deposited 
on solid substrates. 
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